IPJ 2004/01 



Pi of the Sky — all-sky, real-time search for 
fast optical transients*** 

A. Burd a , M. Cwiok b , H. Czyrkowski b , R. Dabrowski b , 
W. Dominik b , M. Grajda a , M. Husejko a , M. Jegier a , 
A. Kalicki a , G. Kasprowicz a , K. Kierzkowski b , K. Kwiecinska c , 
L. Mankiewicz d , K. Nawrocki f , B. Pilecki , L.W. Piotrowski b , 
K. Pozniak a , R. Romaniuk a , R. Salanski a , M. Sokolowski f , 
D. Szczygiel 6 , G. Wrochna f '*, and W. Zabolotny a 

11 Institute of Electronic Systems, Warsaw University of Technology, 
Nowowiejska 15/19, 00-665 Warsaw, Poland 

h Institute of Experimental Physics, Faculty of Physics, Warsaw University, 
Hoza 69, 00-681 Warsaw, Poland 

c Department of Mathematics and Natural Sciences, Cardinal Wyszynski 
University, Dewajtis 5, 01-815 Warsaw, Poland 

d Center for Theoretical Physics, Polish Academy of Science, Al. Lotnikow 32/4-6, 

00-668 Warsaw, Poland 

e Astronomical Observatory, Warsaw University, Al. Ujazdowskie 4, 
00-478 Warsaw, Poland 

{ Soltan Institute for Nuclear Studies, Hoza 69, 00-681 Warsaw, Poland 



Abstract 

An apparatus to search for optical flashes in the sky is described. It has been 
optimized for gamma ray bursts (GRB) optical counterparts. It consists of 2 x 16 
cameras covering all the sky. The sky is monitored continuously and the data are 
analysed on-line. It has self-triggering capability and can react to external triggers 
with negative delay. The prototype with two cameras has been installed at Las 
Campanas (Chile) and is operational from July 2004. The paper presents general 
idea and describes the apparatus in detail. Performance of the prototype is briefly 
reviewed and perspectives for the future are outlined. 
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1 Introduction 



Gamma ray bursts (GRB) are one of the most intriguing phenomena in the 
Universe. Several thousands of them have been observed by gamma ray de- 
tectors carried by satellites. In order to understand their nature, observations 
should be carried out in all wavelengths. Optical domain is especially impor- 
tant, because of high spatial resolution and possibility of detailed spectral 
analysis. Unfortunately, only a few dozens of GRB sources have been seen in 
visible light. Moreover, most of them have been observed as late as hours or 
days after the burst as very faint afterglows. So far, only one GRB was caught 
by an optical device within the first minute after a GRB (Akerlof et al. 1999). 
The observation was made by a small robotic telescope ROTSE equipped with 
photo-lenses of 10 cm aperture (Akerlof et al. 2000). Typical professional tele- 
scopes have too large inertia and too small field of view to be able to react 
promptly to GRB alerts from satellites. 

Since then, many small robotic telescopes have been build and installed around 
the globe in order to search for GRB optical counterparts. Representative ex- 
amples are BOOTES (Ceron et al. 2003), MASTER (Lipunov 2004), RAP- 
TOR (Vestrand et al. 2003), and ROTSE-III (Akerlof et al. 2003). Recent 
review was written by (Andersen & Pedersen 2004). 

Typically, they have optics of rather short focal length and cover relatively 
large field of view (1-100° ). They await for GRB alerts distributed by the 
GCN system (Barthelmy et al. 2000) and move to the target as soon as possible 
in order to take exposures. During long periods between GCN alerts these 
systems typically observe a chosen Field Of View (FOV) and collect data in 
search for transient phenomena of astrophysical interest. 

Unfortunately, so far the GRB930123 recorded by ROTSE remains the only 
optical observation at the time comparable to the duration of the gamma burst 
itself. The reason is two-fold. First, decision making process and propagation 
of the triggering information from a satellite to the observing device takes 
time. Second, the inertia of the devices, although much smaller than that for 
large telescopes, is still responsible for a significant delay. During this time 
the object in question is fading rapidly and it goes beyond the limiting mag- 
nitude of the device before the first exposure begins. It is expected that alert 
system associated with the "Swift" satellite (Gehrels et al. 2004), scheduled 
for launch this year, will operate more efficiently and limit delay related to 
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signal propagation from satellite detector to GCN. However, even the fastest 
systems cannot guarantee systematic observations of the critical region of the 
sky prior to the satellite trigger. Certainly, such observations may reveal addi- 
tional phenomena associated with GRB and shed a new light on their physical 
nature (Paczynski 2001). 



2 The idea of the real-time search with multilevel trigger system 



In order to overcome the two major problems of classical robotic telescopes 
we propose a different approach, based on author's experience from particle 
physics experiments. Trigger propagation time can be eliminated if the device 
has self-triggering capability. Inertia of the system does not matter if the 
object in question is already inside FOV before the investigated phenomenon 
(e.g. GRB) takes place. The apparatus described in this paper exhibits those 
two features. 

We propose to build a system consisting of a number of CCD cameras covering 
as wide field of view as possible. The cameras monitor continuously the sky 
taking relatively short (5-10 s) exposures. The data are analyzed on-line, in 
search for optical transients. The idea is simple: it is enough to check for 
a presence of a star-like object in a given frame, which was not present in 
preceding frames. However, practical realization is difficult, because of large 
data stream involved (see an example below). It is impossible to invent a single 
algorithm, which is fast enough, and has high efficiency and low rate of false 
triggers at the same time. The problem could be solved by implementing a 
multi-level trigger system. The first level algorithms are very simple and have 
high efficiency for interesting events, but they produce a lot of background. 
The rate of background events at this stage can be several orders of magnitude 
higher than the rate of real events we are looking for. The only purpose of the 
first level algorithms is to reduce the data stream to be analyzed by higher 
levels. Thus, the second level algorithm can be somewhat more complicated 
and perform better background rejection. The highest levels deal only with a 
very low rate of suspected events and can employ sophisticated algorithms to 
clean up the final sample. 

Selected events can be submitted to larger telescopes to follow and/or can be 
checked against GRB triggers from other sources, even if they come much later. 
The data cannot be stored for a long time, because of limited disc space, but 
can be temporarily kept in some buffer, to examine any late arriving external 
alerts. 
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3 "7r of the Sky" apparatus design 



Practical realization of the ideas outlined in the previous section may look as 
follows. One observing module consists of a 4 x 4 matrix of CCD cameras. The 
camera is based on CCD with 2000 x 2000 pixels 15 x 15 fim 2 each. Equipped 
with lenses of focal length / = 50 mm, the camera covers 33° x 33° field of 
view. The module of 16 cameras covers almost all visible sky down to 20° 
above the horizon. This is a solid angle above 7r steradians and hence the 
name of the apparatus — "n of the Sk\j\ 

The most common background sources for flash recognition algorithms are 
cosmic rays crossing the CCD and sunlight reflexes from artificial satellites. In 
order to fight this background two twin modules should be installed at some 
distance of the order of 10 km. Cosmic rays can be eliminated by taking simple 
coincidence. Flashing satellites could be resolved by their parallax. 

Gamma ray bursts typically have duration of 0.01-100 s. GRB930123 flash 
observed by ROTSE also had the decay time of the order of dozens of seconds. 
Therefore, the exposure time of 2-10 s seems to be optimal to search for GRB 
related flashes. With the focal length of / = 50 mm and the pixel size of 

15 mm one can take 5 s exposures on a static mount. This is very convenient 
and significantly reduces the cost. The CCD sensitive area is 30 x 30 mm 2 and 
can be illuminated with standard camera lenses, commercially available for a 
low price. 

With 5 s exposures and 1 s readout time one can take 10 frames per minute, i.e. 
5000-10000 frames per night. This imposes severe requirement for the mechan- 
ics of a shutter, which for a few years of operation must sustain 10 7 opening 
cycles. Single frame from one camera occupies 8MB (megabytes). The aver- 
age data rate from a single camera is thus 1.3MB/s = 3GB/h. Two modules, 

16 cameras each, produce the data stream of lOOGB/h, i.e. lTB/night. This 
poses the limit on the temporary data storage. Most of the analysis must be 
done on-line, in real time. Only carefully selected data can be kept for further 
inspection. 

Realisation of the project outlined above is foreseen in the following steps: 

• Phase-0. Tests with one and with two cameras on a fixed mount. 

• Phase- 1. Two cameras observing the same field, installed on a robotic mount 

• Phase-2. Two modules, 16 cameras each, as described above 

• Phase-3. Four pairs of modules installed around the globe, to maximise the 
coverage and observing time 
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4 Phase-0 tests 



The phase-0 of the project has been already completed. Tests have been per- 
formed at Brwinow (52.14725° N, 20.71850° E), 30km West of Warsaw. We 
have started with a modified G-16 camera (Smith 2003) equipped with a 
Kodak CCD KAF-0401E. It has 768 x 512 pixels, 9 x 9//m 2 each. Images 
from about 50 nights between November 2002 and September 2003 have been 
recorded. In the meantime the first prototype of 2000 x 2000 pixel camera has 
been built. It was used for further tests from October 2003. The first final 
camera for the phase-2 has been completed in February 2004, whereas the 
second one in June 2004 (Fig. 1). 

Almost 200 GB of data have been collected in different configurations. Those 
include both, fixed and robotic mounts. The data have been used thoroughly to 
design and tune flash recognition algorithms. Major background sources have 
been identified and their rates have been measured. At least one interesting fast 
optical transient has been observed. The results will be published elsewhere. 




Fig. 1. The "n of the Fig. 2. The Las Campanas site. Left to right: ASAS dome 

Sky" cameras and the housing "n of the Sky" apparatus, ASAS 10" telescope 

mount at the test site dome, Control Room, 
in Poland. 



5 "7r o/ the Sky" phase- 1 system 

The phase- 1 system has been installed at the Las Campanas Observatory 
in Chile in June 2004 (Fig. 2). Regular operation started in July 2004. The 
system consists of two custom designed CCD cameras installed on a robotic 
mount (Fig. 1). 

The cameras are based on CCD442A sensor by Fairchild Imaging. The CCD 
has 2032 x 2032 sensitive pixels, 15 x 15/im 2 each. It is read out with a fre- 
quency 2 MHz/pixel, so the entire matrix is read out in 2 s. After amplification, 
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the signal is digitized by 16-bit ADC and stored in a memory. The camera 
is read out and controlled by a PC through a fast USB 2.0 interface. Data 
transfer takes less than one second and can be done while the next exposure 
is already being taken. Readout speed, amplifier gain and other parameters 
are programmable via USB. The sensor is cooled with a stack of two Peltier 
modules about 35 degrees below the ambient temperature. Special heavy-duty 
mechanical shutter was designed to sustain over 10 7 opening cycles. A proto- 
type has been tested at high frequency in a lab. The first signs of degradation 
have appeared after 1.2 • 10 7 cycles. Details of the camera design will be pub- 
lished elsewhere. 

Each camera is equipped with Planar-T* photo-lenses by Carl-Zeiss of a focal 
length / = 50 mm and an aperture d = //1.4 . Focusing is performed by a 
step motor with a controller build into the camera and remotely controlled 
through the USB. The effective field of view is 33° x 33°. The two cameras 
are installed on a common mount in such a way that they observe the same 
field. The mount we use was originally designed for the ASAS experiment 
(Pojmanski 2000) and modified to suit the "n of the Sktf needs. 

It is driven by step motors controlled by a PC through RS232 serial interface. 
The mount can reach any point in the sky in less than one-minute time. 
The whole setup is installed in the ASAS dome together with other ASAS 
telescopes. 



6 Operation and data flow 

The apparatus is controlled by a PC located inside the dome. Second PC, 
located in a nearby Control Room is used for off-line data analysis. The system 
is fully autonomous, but also fully controllable via Internet. During the normal 
operation the system runs autonomously according to the preprogrammed 
schedule. Dedicated script language has been developed to make the schedule 
programming easy and flexible. 

For most of the time the cameras follow the field of view of the HETE satel- 
lite (Ricker et al. 2002). Its position is read out from the Internet in regular 
intervals and the mount position is automatically corrected accordingly. If the 
HETE FOV is not visible, another location in the sky is programmed. The 
system is also listening to GCN alerts received directly and through another 
server in Warsaw as a backup. Should an alert located outside the current FOV 
arrive, the mount automatically moves towards the target and exposures are 
being taken. Twice a night an all sky scanning is performed. 16 fields are vis- 
ited and three images of each are taken by both cameras. A single scan last 
about 20 minutes. 
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Because the mount follows the sideral movement, one can take longer expo- 
sures than in the case of a fixed system. We have chosen 10 s as a compromise 
between the magnitudo reach and time resolution for short flashes. The im- 
ages are immediately analyzed while in the computer RAM in search for flashes 
with a rise time of the order of seconds. Then, they are temporarily stored on 
a disc and can be reexamined in case of late arrival of an external alert. If a 
flash candidate is found the 100 x 100 pixel samples of ±7 frames are stored 
permanently for the record. 

In the meantime, the images are copied to the second PC, which superposes 
the images and searches for optical transients with a rise time of minutes. 
During the day, two analyses are performed in parallel on the temporarily 
stored data. The first PC runs fast photometry on individual frames, which 
can be used later to study rapidly varying objects. The second PC performs 
precise photometry (Pojmanski 2000) on images superposed by 20. This could 
be used to study variable stars etc. The results are stored permanently on a 
disc. Out of almost 30 GB of data taken every night, about 2 GB of results 
is stored permanently. After 2-3 months a 200 GB removable disc with the 
results is replaced and taken to Warsaw for further analysis. 



7 Software 

Both PC are running the Linux operating system. The software is mostly 
custom written in C++. It consists of a number of modules taking care of 
different devices: the mount, the cameras, the data acquisition system, the 
GCN server, etc. All modules are governed by the central module "piman" 
scheduling the tasks and controlling the information flow. The communica- 
tion between modules it based on C0RBA. Image processing classes make use of 
"cfitsio" package (Pence 1993). Flash recognition algorithms and fast pho- 
tometry are custom developed, whereas precise photometry and astrometry is 
adopted from ASAS (Pojmanski 2000). 

Monte Carlo technique was extensively used in development of the algorithms. 
Optical flashes were simulated by cutting star images from real sky images and 
pasting them into images under study at random locations. This method was 
used to evaluate efficiency of flash detection algorithm. Details are described 
in (Piotrowski & Sokolowski 2004). 

The flash recognition algorithm compares a given frame with several preceding 
ones. It searches for a star-like object, which is missing on preceding frames. 
Most of the false events are caused by cosmic rays, but those are easily elimi- 
nated by coincidence of the two cameras. The most severe background is due 
to the sunlight reflexes from artificial satellites. We try to eliminate those in 
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two ways. First, we search for aligned flashes in a single frame, or in different 
frames. Second, we check the time and location of the flash against a database. 
Every evening a fresh database is built by merging several ones available on 
the Internet (McCants 2004). 



8 System reliability 

The apparatus operates at the Las Campanas Observatory without a per- 
manent human supervision. Therefore, the system must be very reliable. It 
was achieved by employing self-diagnostics and remote monitoring, as well as 
hardware redundancy and flexible configuration. 

Special care has been taken to ensure seamless recovery after a system failure 
without human intervention at the site. The cameras have a hardware watch- 
dog build in, which automatically resets the camera and thus reestablishes 
connection with the PC in case of protocol failure. In addition to the direct 
Internet connection the two PC are connected together with a Gbit Ethernet. 
It ensures fast data transfer between the two PC and also serves as a backup in 
case one of the direct connections breaks. Both PC have "Wake on LAN" and 
"Boot from LAN" capabilities and can be started from the net in case of a file 
system failure. Remote file system could be installed and used to recover the 
local one. Each computer can be remotely reset or powered down/up by the 
relays driven by the second one. This feature has proved itself to be extremely 
useful in case of system hang-ups. 

The system is designed to be immune to network failures. It can run au- 
tonomously for many days. Each night it reads the current schedule from a 
script and a default script is executed if no current script is found. The system 
can be effectively monitored with a very small bandwidth. Direct communi- 
cation with the system is provided by a console module "pishell" which can 
run at any place and interact with the system exchanging short COBRA packets. 
Basic information about the system (< 20 kB) is automatically copied every 
10 min to a WWW server in Warsaw. Selected sky images are compressed to 
4kB JPEG and also copied to the Warsaw server every 20 min. Should any 
failure occur, the system sends an SMS with the appropriate information to a 
mobile phone of a person in charge. 



9 System performance 

The phase-1 "n of the Skif system operates regularly from July 2004. Till 
the end of October, three nights were lost for the maintenance after a double 
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disk failure and several more due to weather conditions. About half a million 
sky images have been taken, accounting for 4 terabytes of data. Each frame 
contains about 20 000 stars and the results consists of 10 billion photometric 
measurements. 



The limiting magnitude is 10 — ll m for single frames and 11 — 12 m for frames 
coadded by 20 (Fig. 3 and 4). Exact limit depends on the sky background, 
which strongly varies with the Moon phase. It also depends on the position in 
the frame. The corners are somewhat less illuminated than the center. 
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Fig. 3. Fast photometry magnitudo of 
stars vs position in the frame (degrees) . 



Fig. 4. Precise photometry magnitudo 
of stars vs position in the frame (de- 
grees) . 



During the first three months of the regular operation, several short optical 
flashes have been detected. Most probably they are caused by sun reflexes from 
artificial satellites, although they do not correspond to any satellite listed 
in commonly available databases. One optical transient lasting > 10 s was 
observed during tests with a single camera, which cannot be explained by 
a satellite reflex. So far, no optical transient coinciding with a known GRB 
was found. Limits for optical counterparts have been given for a few GRB 
(Cwiok 2004). Detailed results will be published elsewhere. 



10 Perspectives 



Current efforts in the project concentrate in two areas. First, the software 
is still being developed to improve flash recognition algorithms and increase 
precision of the photometry (Fig. 5 and 6). The work is going on to develop a 
database with photometric measurements, which can be used to study variable 
stars, etc. Second, the design of the phase-2 apparatus (2 x 16 cameras covering 
all the sky) is moving from a conceptual to a technical level. At the same time 
we search for optimal locations of the phase-3 systems. 

The "n of the Ski/' apparatus can be also considered as an Earth based pro- 
totype of a lighter, yet more powerful system, which can operate in the space 
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FAST PHOTOMETRY 



PRECISE PHOTOMETRY 





Fig. 5. Precision vs brightness for the Fig. 6. Precision vs brightness for the 
fast photometry. precise photometry. 

(Tsarevsky 2002). On the other hand it stands as a milestone towards future 
telescope "farms" monitoring all sky in search for rare and rapid phenomena. 
The system is already integrated with GCN (Barthelmy et al. 2000) and OTC 
(Kuvshinov & Lipunov 2004) networks. "7r of the Shtf triggers will be soon 
automaticaly sent to the robotic telescopes of ASAS. We also consider to in- 
tegrate our system with more sophisticated networks, e.g. TALON (White et 
al. 2004). 

Technological progress in image sensors, electronics and computers will un- 
doubtly push observational astronomy towards huge data streams, which have 
to be analyzed on-line. Real time event selection mechanisms, large system is- 
sues, maintenance and operation problems are the areas that need a lot of 
research and development. The "ir of the Sktf experiment is one of the first 
steps on this new land. 
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